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ABSTRACT 

We study the galactic population of double compact objects (NS-NS, BH-NS, BH-BH binaries) to 
investigate the number (if any) of these systems that can potentially be detected with LISA at low 
gravitational-wave frequencies. We calculate the Galactic numbers and physical properties of these 
binaries and show their relative contribution from the disk, bulge and halo. It is found that, although 
the formation efficiency of double compact objects is very small, a few thousand of such systems should 
be detectable with LISA depending on the instrument low-frequency sensitivity and the efficiency of 
data analysis techniques for identifying highly eccentric systems. The vast majority of these binaries will 
add up to foreground noise arising from the galactic double white dwarf binaries, and only a very few 
1 — 10) are predicted to stand out above the noise and be resolvable. Detailed population synthesis, 
combined with an updated Galactic model and a realistic LISA signal simulator are used to obtain the 
results. We outline the prospects of the science that could be done with future LISA observations of 
double compact objects. In particular, we note that LISA observations could be used to determine the 
fate of massive binaries that are the most likely progenitors of BH-BH systems and shed some light on the 
common envelope evolution involved in the formation of double compact objects. Finally, the galactic 
double compact object populations presented here may readily serve in LISA data mock challenges^ 



''The data can be obtained upon request from Chris Belczynski. 
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1. INTRODUCTION 

The Laser Interferometer Space Antenna (LISA) is a 
next generation instrument to search for and observe grav- 
itational radiation (e.g., Hughes 2006). A system of three 
satellites (5 million kilometers apart), orbiting the Sun, 
will form an interferometer sensitive to low frequency grav- 
itational radiation (GR) (-- 5 x 10"^ - 1 Hz). The mam 
sources of GR at these low frequencies are inspirals of 
super-massive black holes in the center of merging galaxies, 
extreme mass ratio inspirals of stars and stellar mass black 
holes into the super-massive black holes, and also nearby 
(mostly Galactic) compact binaries. Double white dwarf 
binaries (WD- WD) are the largest population of Galactic 
systems that are expected to produce a confusion-limited 
noise in the detector, with several thousand of the louder 
systems being potentially resolved. A number of stud- 
ies have concentrated on studies of double white dwarfs 
in the context of low- frequency LISA observations (e.g., 
Hils, Bender & Webbink 1990; Farmer & Phinney 2003; 
Nelemans, Yungelson & Portegies Zwart 2004; Ruiter et 
al. 2009a, b). Here, we focus on the other, and much less 
studied group of compact systems; double neutron star 
binaries (NS-NS), black hole neutron stars (BH-NS), and 
double black holes (BH-BH). Although much less com- 
mon than double white dwarfs, these systems may produce 
much stronger GR signals, since these compact objects are 
much denser and more massive than white dwarfs. 

Out of all double compact objects only NS-NS have 
so far been observed. A handful of NS-NS were discov- 



ered in radio surveys (e.g., Lorimer 2005). Seven systems 
are found in the Galactic field, while one belongs to the 
Galactic globular cluster M15 (PSR B2127-11C). Out of 
the seven, three are wide systems with merger time over 
Hubble time (PSR J1518-l^4804, PSR J1811-1736, PSR 
J1829+2456). The remaining four constitute the close NS- 
NS population with relatively short (^ 0.1 — 3 Gyr) merger 
times (PSR B1913-I-16, PSR B1534-M2, PSR J1756-2251, 
PSR J0737-3039). Also one additional field system (PSR 
J19G6-I-0746) is considered as a very probable close NS-NS 
(merger time ~ 0.3 Gyr) binary (Kasian 2008), although 
the possibility that this is a NS-WD binary can not yet be 
excluded. 

Mergers of all types of double compact objects 
are expected to be the prime candidates for high- 
frequency ground based gravitational- wave interferometers 
like LIGO or VIRGO, while mergers of NS-NS, BH-NS and 
BH-NS were proposed as potential progenitors of short- 
hard Gamma-ray bursts (Paczynski 1986). Double com- 
pact object binaries were studied extensively over last 2 
decades, bringing new developments and understanding of 
their formation in particular in context of population syn- 
thesis studies (e.g., Tutukov & Yungelson 1994; Lipunov, 
Postnov & Prokhorov 1997; Portegies Zwart & Yungel- 
son 1998; Bethe & Brown 1998; Belczynski & Bulik 1999; 
Fryer, Woosley & Hartmann 1999; Bloom, Sigurdsson & 
Pols 1999; Nelemans, Yungelson & Portegies Zwart 2001; 
Hurley, Tout & Pols 2002; Voss & Tauris 2003; Dewi 
& Pols 2003; Pfahl, Podsiadlowski & Rappaport 2005). 
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Merger rates were recently presented and discussed by 
Kim, Kalogera & Lorimer (2006; empirical estimates) and 
Belczynski et al. (2007; population synthesis). 

In this study, we analyze the entire Galactic popula- 
tion (disk, bulge and halo) of double compact objects. 
Using population synthesis methods we predict numbers 
and physical properties of double compact objects, then 
wc calculate their spatial distribution and finally estimate 
the low frequency gravitational radiation signal that may 
arise from these binaries. In § 2 we describe our method of 
calculation, in § 3 we present results, in § 4 we compare our 
findings with previous work, and finally in § 5 we present 
our conclusions. 

2. MODELING 

2.1. Population Synthesis 

We have used the population synthesis code StarTrack 
to calculate the numbers and properties of double compact 
objects. The full description of the code can be found in 
Belczynski, Kalogera & Bulik (2002) and Belczynski ct al. 
(2008). The code utilizes a set of stellar models (Hurley 
et al. 2000) that allow for evolution of stars at differ- 
ent metallicities. The compact object formation follows 
self-consistently from the stellar models, extended to the 
formation of the FeNi core (Timmcs, Wooslcy & Weaver 
1996). During the core collapse the fall back and direct 
black hole formation is accounted for (Fryer & Kalogera 
2001) and the newly born compact objects receive natal 
kicks (Hobbs et al. 2005). Formation of low mass neu- 
tron stars throiigh electron captiirc siipcrnovae is also ac- 
counted for (e.g., Podsiadlowski et al. 2004). Binary in- 
teractions are treated in detail and the various processes 
were calibrated using either results of detailed evolution- 
ary calculations (e.g., Wellstein, Langer & Braun 2001 for 
mass transfer sequences) or specific sets of observations 
(e.g., Levine, Rappaport & Zojcheski 2000 for tidal inter- 
actions). 

All of our popiilation synthesis calculations implement 
the standard evolutionary model presented in Belczynski 
et al. (2008). We employ our standard model to evolve 
three different Galactic populations: disk, bulge and halo. 
The dynamical evolution of stars in these components is 
described in the following section. Here we describe the 
evolutionary aspect of the stars in each component. The 
disk is assumed to have a stellar mass of 4.0 x 10^° M©, 
stars have solar metallicity Z = 0.02 and we assume the 
age of the population to be 10 Gyr with constant star 
formation throughout the disk lifcitimc (i.c;.. star forma- 
tion rate is 4 Mq yr^^). The bulge is assumed to have a 
stellar mass of 1.1 x 10^° Mq, stars also have solar metal- 
licity Z = 0.02 and we assume the age of the popula- 
tion to be 10 Gyr with a burst of star formation that 
lasted throughout first Gyr of bulge evolution (i.e., all the 
stars were born between 9 and 10 Gyr ago; a uniform dis- 
tribution of birth times was adopted for all the stars in 
this time interval). The halo is assumed to have a stellar 
mass of 0.1 X 10^° M©, the stars have sub-solar metallicity 
Z = 0.001 and we assume the age of the population to be 
13 Gyr with a burst star formation at the very beginning 
of halo evolution (i.e., all the stars were born 13 Gyr ago 
in a delta function-like burst). 

We assume that in each population the binary fraction is 



50% (2/3 of stars are in binaries, and 1/3 are single stars). 
We do not consider any dynamical interactions between 
stars; i.e., evolution of stars in globular clusters is not 
accounted for and we evolve only field populations. How- 
ever, we note that despite their relative small stellar mass 
content globular clusters may contribute significantly to 
formation of BH-BH binaries (Kulkarni, Hut & McMillan 
1993; Portcgics Zwart & McMillan 2000; Guhekin, Miller 
& Hamilton 2004; O'Leary et al. 2006; Sadowski et al. 
2008). On the other hand, formation of double compact 
object systems with neutron stars was found to be rather 
inefficient in globular cluster environments (e.g., Phin- 
ney 1991; Grindlay, Portcgics Zwart, & McMillan 2006; 
Ivanova et al. 2008). The initial mass of single stars 
and binary primaries (more massive component) are drawn 
from the three component broken power-law Initial Mass 
Function (IMF) . We employ the IMF from Kroupa, Tout & 
Gilmore (1993) and Kroupa & Weidner (2003) with slope 
of -1.3 (0.08 < M < 0.5 M0), -2.2 (0.5 < M < 1.0 Mq), 
-2.7 (1.0 < M < 150 Mq). The binary secondaries are 
drawn through flat mass ratio distribution within the en- 
tire mass range (0.08 < M < 150 M©). We follow only the 
evolution of massive binaries as only these can form dou- 
ble compact object binaries. We require that the mass of 
the primary and secondary are the following: Mi > 6 Mq 
and M2 > 4 Mq. Even though it is expected that only 
stars more massive than ~ 8 Mq can form compact objects 
(i.e., neutron stars and black holes), the binary evolution 
(through mass transfer and rejuvenation) can alter such a 
simple picture. The binary initial separations are drawn 
from distribution that is flat in logarithm (Abt 1983) with 
maximum separations reaching 10^ R©. Binaries are al- 
lowed to be eccentric, and we use thermal distribution for 
initial eccentricities (i.e., 2e; Heggie 1975; Duquennoy & 
Mayor 1991). 

It was noted by Belczynski et al. (2007) that many pro- 
genitors of double compact object binaries evolve through 
a common envelope phase while a donor star is crossing 
Hertzsprung gap. Such a star does not have a well de- 
veloped core-envelope structure (e.g., Taam & Sandquist 
2000) and once inspiral in a common envelope is started it 
may very likely lead to a merger whether there is enough 
orbital energy to expel the envelope or not. If such merg- 
ers are taken into account they lead to a significant de- 
crease in the formation of close double compact objects. 
To estimate the impact of this uncertainty on the LISA 
signal we perform two calculations. In one (Model A) we 
allow for survival in the case of a common envelope with 
a Hertzsprung gap donor (i.e., standard energy balance is 
tested to check for system survival; Webbink 1984), while 
in another (model B) we assume that all such common 
envelope events lead to a merger aborting further binary 
evolution and thus depleting the population of double com- 
pact objects. 

2.2. Galactic Model 

The binaries are distributed at birth throughout the 
Galaxy according to stellar density models. We assume 
that the densities are independent of time. For the bulge, 
we choose a spherical density with a normal distribution 
in the radial coordinate and a cut-off radius of 3.5 kpc, so 

Pb oc e-('-/'-°^)', (1) 
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where rpb = 500 pc (Nelemans et al. 2004). The disk pop- 
ulation is assumed to be axially symmetric with cyhndrical 
radius and vertical height distribution given by a double 
exponential: 

Pd oc e-■"■/■""^e-l^l/"''^ (2) 

with i?od = 2.5 kpc and zod — 200 pc. Finally, we dis- 
tribute the halo binaries according to a simplistic spherical 
model with: 



1 + ^ 



-3.5 



(3) 



with oo = 3.5 kpc. The disk and bulge systems are given 
initial rotational velocities in the plane of the disk, and 
the halo systems are assumed to have circular orbits corre- 
sponding to their initial positions. All systems propagate 
through a Galactic potential and their trajectories change 
due to the kicks received at the birth of each compact 
object. 

The Galactic gravitational potential is the sum of the 
bulge, disk and halo potentials. The disk and bulge are 
described the Miyamoto and Nagai (1975) potentials: 



$(i?,z) 



GM, 
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where R — \f x" 



r.2 -i-y2^ |-]-^g jjif^gx i corresponds to bulge 
and disk, and hi are the parameters, and Mi is the 
mass. The halo potential is assumed to be 
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where where Tc is the core radius, and Mh is the parame- 
ter describing the mass of the halo (Paczynski 1990). The 
mass of such halo is divergent so we introduce a cutoff ra- 
dius Tcut beyond which the halo density falls to zero, and 
the potential is $(r) cx r^^. The bulge is described by 
Ml = 1.12 X 10^° Mq, ai = Okpc, and hi = 0.277kpc, 
and the disk by M2 = 8.78 x 10^° Mq, 02 = 4.2 kpc, and 
62 — 0.198 kpc. For the halo potential we use = 6.0 kpc, 
rc„t = 100 kpc, and Mu = 5.0 x 10^° Mq. We note that 
the masses of the potentials differ from the masses of the 
stellar components in Section [2.11 because we include non- 
stellar matter and dark matter in the potentials. 

2.3. LISA Signal Simulator 

We calculate the expected signal in LISA as the Michel- 
son signal from an equal-arm interferometer using the 
long wavelength approximation for frequencies below 3 
mHz (Cutler 1998; Benacquista, DeGoes & Lunder 2004) 
and the rigid adiabatic approximation (Rubbo, Cornish 
& Poujade 2004) above this frequency. The gravitational 
waveform for an eccentric binary with angular frequency 
bj = 27r/Porb is calculated in the quadrupole approxima- 
tion (Peters & Mathews 1963; Peters 1964). We use the 
form given in Pierro et al. (2001): 

h+ — —j=\Jy\ + cos^ -d) \2,hxySm2Lp + hx-yC0s2Lp\ 
v2 

- (1 - COS^ §) h, + y\ (6) 

COS "3 

hx — — [2hxy cos 2(p — h^_y siii2(p] (7) 
v2 



where •& and if are the usual spherical angular coordinates 
in the frame centered on the binary with the x-axis along 
the periastron. The metric components appear as: 



h^"^ sin {nujt + n(j)o) 



71=1 



hx±v = Ir^ly cos (nujt + ncpo) 



(8) 
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where (jjQ is the phase of the binary at i = and 



"'xy 



An) 



(n) 

x+y 



ho n\/l - e2[J„_2(ne) + Jn+2{ne) 
-2J„(ne)] 

2ho {n [Jn-2{ne) - Jn+2{ne)] 
-2ne [J„_i(ne) - J„+i(ne)] + 2J„(ne)](ll) 
-4/ioJ„(ne) 



(10) 
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where Jn{x) are Bessel functions. The overall amplitude 
of the gravitational wave is proportional to: 



ho = 



2G5/3,2/3^,5/3^ 



(13) 



where r is the luminosity distance from the earth to the 
binary. We assign each binary an arbitrary orientation, 
and so the contribution of each binary to the LISA sig- 
nal is described by 11 parameters. The masses (Mi and 
M2), luminosity distance (r), and orbital frequency (/orb) 
are used to construct the overall amplitude and the rel- 
ative amplitudes of the harmonics are obtained from the 
eccentricity (e). The sky location is obtained from the fi- 
nal positions {9s and (f>s) of the binaries after propagation 
through the Galactic potential. The arbitrary orientation 
includes the initial phase {(f>o), the argument of the perias- 
tron (7), and the direction angles of the angular momen- 
tum vector {01 and (j>i). We compute the one- year time 
domain signal for each harmonic in the waveform up to 
a frequency of 30 mHz. We have not included periastron 
precession, although this may prove to be important for 
longer observation times. 

For e = 0, all the power in the waveform is concentrated 
at the n = 2 harmonic. At non-zero eccentricity, the power 
in this waveform is spread out over several harmonics of 
the orbital frequency. If we define nmax(e) to be the har- 
monic with maximum power for eccentricity e, then for 
e ^ 0.8, nmax(e) can be approximated by 



1.6 



(l-e)i-5' 



(14) 



The factor of 1.6 arises from maximizing g{n, e) from Pe- 
ters & Mathews (1963), but can be understood by realizing 
that the speeds of the components are larger at pericen- 
ter than for the corresponding circular orbit. This can be 
understood by taking into account that the gravitational 
waves must emit close to the frequency of a circular orbit 
with radius equal to the periastron, rperi, which is given 
by: 



friiax C>C ^pcri (1 



e) 



/orb ^ ^max/orb- 



(15) 
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Therefore, binaries with /orb below the low-frequency limit 
(/crit) for LISA sensitivity may still be observable by LISA 
if '^max/orb > /crit- We use the following approximate cri- 
terion to separate potential LISA sources from the entire 
double compact population 

Porb < 1.6(1 - e)-i Vcri't (16) 

where /crit = 5 x 10^^ Hz. 

The LISA noise is simulated by assuming the power 
spectral density of the noise is made up of position (or 
shot) noise given by (Cornish, 2001): 

S-np = 4.8 X 10"''2 Hz-\ (17) 

and an acceleration noise (converted to strain) given by: 

= 2.3 X 10-^0 (^^7^)' Hz-^ (18) 
These separate components are combined according to: 

Sn = 45r.p + 85„a (l + COS^ (///,)) , (19) 

where /* = c/2ttL with the armlength of LISA taken to 
be L = 5 X 10^ m. We roll off the acceleration below /crit, 
so that S'„a(/ < /crit) = 5'na(/crit)- In reality, the LISA 
noise will probably not follow this simple power law all the 
way down to our choice of /crit , but will begin to rise at a 
higher frequency. 

3. RESULTS 
3.1. Physical Properties & Numbers 
3.1.1. Numbers: Model A 

In Table 1 we present the number of double compact 
objects predicted for the present time in our Galaxy. The 
numbers are given for each Galactic component (i.e., disk, 
bulge, halo) and the contribution of each type of double 
compact object binary (i.e, NS-NS, BH-NS, BH-BH) is 
shown. Additionally to the total number of Galactic dou- 
ble compact objects, we list the number of LISA systems. 
The LISA systems are the binaries that are most likely to 
have right properties (see eq. [TB]) to produce a GR signal 
in the frequency sensitivity band of LISA. In the following 
we discuss only results for Model A. 

As a whole, the Galactic population of double compact 
object is dominated by disk (66%) and bulge (30%) sys- 
tems with a small contribution of halo binaries (4%). This 
reflects the fact the the number of double compact objects 
is proportional (for a given evolutionary model) to the stel- 
lar mass. Since the disk stellar mass is the largest, the disk 
population of double compact objects is the most numer- 
ous. The bulge stellar mass is about half that of the disk, 
and the number of bulge double compact object binaries 
is about half of those in the disk. Similarly for the halo. 
The numbers are not exactly proportional to the initial 
stellar mass, since the star formation history and initial 
conditions (e.g., metallicity) alter the formation efficiency. 
However, these factors are only secondary in comparison 
to the dependence of double compact object number on 
the initial stellar mass. 



Within the entire double compact object population, 
BH-BH systems are dominant (67%), with a significant 
contribution of NS-NS binaries (28%) and the very small 
relative number of BH-NS systems (5%). The dominance 
of BH-BH systems is a consequence of the natal kick model 
employed, i.e., neutron stars are given full kicks, while 
black holes (due to formation through either partial fall 
back or direct stellar collapse) receive smaller or no natal 
kicks at all. This naturally leads to the survival of BH-BH 
progenitors (especially those on wide orbits) while the vast 
majority of NS-NS progenitors are disrupted at either the 
first or second supernova. Progenitors of BH-NS systems 
are depleted by two processes. First, the large difference 
in initial mass of two stars (one very massive to make a 
BH, and the other of lower mass to produce a NS) tends 
to lead to a merger in the first RLOF episode early in 
the evolutionary history of such systems. Second, many 
systems that survive the first RLOF and subsequent for- 
mation of a BH, are disrupted by the natal kick during NS 
formation. 

The population of potential LISA sources is quite dif- 
ferent than the entire double compact object population. 
These are the binaries that were formed at close orbits 
(^orb ^150 day), and they constitute only a smaU frac- 
tion of the entire population. The NS-NS LISA binaries 
are 30% of entire Galactic NS-NS population, while it is 
7% for BH-NS and only 1% for BH-BH LISA binaries. The 
decreasing contribution to the LISA population is due to 
the fact that many more BH-BH systems were allowed to 
form on rather wide orbits while wide NS-NS binaries were 
mostly disrupted by natal kicks, as explained above. 

The bias towards disk binaries is even more pronounced 
for the LISA binaries, in that they mostly originate from 
the disk (78.8%), with significant contribution from the 
bulge (20.9%) and a negligible number of the systems in 
halo (0.3%). This increased disk contribution (relative to 
the entire double compact object population) is an effect 
of star formation history. Basically, the LISA group con- 
sists of such close binaries that they can merge (and disap- 
pear from the population) in less than the age of Galaxy. 
Since the disk is forming stars in an approximately con- 
stant manner there is a constant supply of close binaries, 
while for bulge and halo that have formed all the stars 
~ 10 Gyr ago, many close binaries have already merged. 

The most striking feature of LISA population is that it 
mostly consists of NS-NS binaries (90.4%), with only mi- 
nor contributions from BH-NS (3.7%) and BH-BH (5.9%) 
systems. For single star evolution with the adopted IMF 
it is predicted that NSs outnumber BHs by about factor 
of 5. 

Binary evolution further affects the numbers. Double 
neutron star progenitors are more likely to avoid merg- 
ers in RLOF episodes as the two stars in the binary are 
close in mass, while for BH-BH progenitors as mentioned 
above it is likely that dynamical instability develops and 
a progenitor system enters a common envelope that may 
lead to a merger. Since, in the case of the LISA pop- 
ulation, we consider only the population of progenitors 
that start in close orbits (to form the very close double 
compact objects) kicks that are higher for neutron stars 



''The IMF slope is -2.7, neutron star progenitor mass range in our single star model is 7.5 — 21 Mq while for black holes it is 21 — 150 Mq 
(see Belczynski et al. 2008; their Fig. 1). 
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are not an important factor that could allow binaries with 
BHs to outnumber NS-NS systems. Additionally, neutron 
stars with less massive (and smaller in size) progenitors 
can start their evolution from closer orbits and therefore 
form double compact object binaries on closer orbits as 
compared with BH progenitors. The dominance of NS- 
NS binaries over BH-NS and BH-BH systems was already 
noted in our models for a different Galactic population (bi- 
naries with even closer orbits) of double compact object 
mergers (Belczynski et al. 2007). 

3.1.2. Physical Properties: Model A 

In Figure[l]we show the characteristic properties of dou- 
ble compact objects that may be detected by LISA. The 
panels show distributions of orbital frequencies (/orb = 
1/Porb; top), eccentricities (middle) and chirp masses (bot- 
tom) of LISA double compact objects. There is a subtle 
difference between the distributions presented here and in 
various papers for double compact objects (e.g., Belczyn- 
ski et al. 2007). Mostly the results for double compact 
objects are calculated for a uniform population of stars 
and are presented at the double compact object formation 
or for the population that will merge within some spe- 
cific time (e.g., Hubble time). The distributions presented 
here, however similar they may appear, are none of the 
above. First, our distributions originate from the mixture 
of stellar populations used to represent the Galaxy (dif- 
ferent metallicities and star formation histories). Second, 
double compact objects form at a given time (set primarily 
by a specific star formation history) and then are evolved 
to the present time. During the evolution each compact 
object is subject to the emission of GR and its orbit (pe- 
riod and eccentricity) is decaying. Finally, we only select 
and present double compact objects that satisfy Eq. [121 

The frequency distributions are similar for all types 
of double compact objects and they span a wide range: 
/orb - 10"^ - 10"^ Hz (Porb 0.02 - 150 day) and peak 
at /orb 5 X lO^'^ Hz (Porb ^ 0.5 day). The shape of the 
frequency distributions may be understood in the follow- 
ing way. For frequencies lower than /crit = 5 x 10~^ Hz (or 
large orbital periods) the number of systems is falling off 
since the eccentricity of a given system needs to increase 
with decreasing frequency in order for the system to be de- 
tectable (see eq.fTB)). However, the number of systems with 
higher and higher eccentricities is gradually decreasing de- 
spite the fact that eccentricity distributions are rather fiat 
(see below). For high frequencies (/orb > /crit) the systems 
are so tight (Porb < 0.5 day for a circular binary) that the 
orbital decay is very fast due to a very effective emission 
of GR. The tighter the system (higher /orb), the shorter 
time it takes for a system to merge causing a depletion in 
number with an increasing GR frequency. 

The eccentricity distributions are rather flat for all dou- 
ble compact objects. For NS-NS and BH-NS binaries the 
distributions are slightly skewed toward the high e-values, 
while it is just the opposite for BH-BH systems. The ec- 
centricity distribution is the direct result of the second 
supernova asymmetry; both natal kick and mass loss from 
the exploding star. Since systems with neutron stars re- 
ceive, on average, larger second kick the NS-NS and BH- 
NS binaries are more eccentric than BH-BH systems. The 
rather high fraction (~ 40%) of BH-BH systems with small 
eccentricities (e < 0.2) is either the result of direct collapse 



of a star to a BH or significant fall back of material during 
the second supernova explosion. 

Chirp mass (Mchirp = (MiJVhf-^/iMi + Ma)^-^) distri- 
butions are very different for the three subclasses of dou- 
ble compact objects. The NS-NS binaries have a chirp 
mass distribution that peaks at Mchirp ~ 1-2 Mq with 
a tail that extends to Mchirp ~ 2 Mq. Systems with 
black holes have much fiatter distributions but spanning 
a wide range of chirp masses; Af chirp ^ 1.5 — 4 Mq for 
BH-NS and Mchirp ^ 2.5-9 Mq for BH-BH binaries. 
Most neutrons stars form with mass 1.2 — 1.4 Mq either 
through electron capture supernova, or through regular 
FeNi core collapse (stars with initial mass Mzams < 18 Mq; 
Timmes et al. 1996) and they explain the peak at low 
chirp mass. The accretion onto the first born neutron 
star was calculated (Belczynski et al. 2007) to be rather 
small (~ 0.05 — 0.1 Mq) and just high enough to mildly 
recycle the pulsar. The mild recycling and the range 
of neutron star masses predicted here agrees rather well 
with the observed masses in the known NS-NS systems, 
all of which are found with 1 — 1.5 Mq (e.g., Lorimer 
2005). Only very few neutron stars are predicted with 
higher masses in the NS-NS population since they are 
formed from initially massive stars and in a narrow range 
(18 < Mzams <^ 21 M0) and so they are not predicted to 
show up in such a small observed sample (8 systems) . For 
BH-NS and BH-BH systems the distributions are mostly 
shaped by the mass distribution of black holes. The black 
holes are predicted within the very wide mass range and 
that makes the chirp mass distributions wide as well, wider 
for BH-BH systems and somewhat narrower for BH-NS 
binaries since for these systems one component is always 
of a low mass (^ 1.3 Mq). The lowest mass black holes 
are found just over 2.5 Mq as we adopted such a value 
as the maximum neutron star mass and they can be as 
massive as ~ 15 Mq for disk and bulge (high metallic- 
ity) while they can reach even higher masses ~ 30 M0 for 
halo population (low metallicity environment). The very 
high masses of black holes predicted for the halo are the 
result of the evolution of stars in a low metallicity environ- 
ment; stars lose much less mass in winds so there is a much 
larger mass reservoir for compact object formation as com- 
pared to high metallicity models (disk and bulge). Here 
our predictions are also consistent with the observations of 
binaries hosting confirmed stellar mass black holes. Stel- 
lar mass black holes are found in our Galaxy with masses 
reaching ~ 15 Mq (e.g., Casares 2006), while in the ex- 
ternal galaxy IC-10 with a low metallicity environment a 
black hole was found with mass ~ 24 — 33 Mq (Prestwich 
et al. 2007; Silverman & Fflipenko 2008). Note that al- 
though we predict such high mass black holes in the halo 
of our Galaxy there are very few of them, and they do not 
have very massive companions as the highest chirp masses 
for BH-BH binaries are Mchirp < 10 Mq. The system in 
IC-10 was predicted to form a BH-BH binary with chirp 
mass of A/chirp ^ 20 Mq but with a rather short coales- 
cence time ~ 2 — 3 Gyr (Bulik, Belczynski & Prestwich 
2008). Therefore, if such a system has formed in the halo 
of our Galaxy it has most probably merged by now. 

3.1.3. Model B 
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In Table 1 we also list the numbers of double compact 
object binaries for our Model B, in which we do not allow 
for survival through common envelope with Hertzsprung 
gap donors. 

For entire population the most affected systems are BH- 
NS systems (reduction by factor of ~ 2), and BH-BH bina- 
ries (reduction by ^ 1.4) while NS-NS binaries are affected 
the least (reduction by ~ 1.2). All these changes are not 
very large especially in the light of the limited prospects 
for detecting the majority of these binaries; since the ma- 
jority of these systems are too wide to ever make it to the 
LISA frequency range or to merge within a Hubble time 
to be detected by ground-based detectors like LIGO or 
VIRGO. 

For LISA binaries, the differences are much more pro- 
nounced. BH-BH binaries are reduced by factor ~ 250, 
BH-NS by factor ~ 10, while NS-NS by factor - 1.5. The 
progenitors of close LISA binaries are subject to one or 
more common envelope events and thus they are greatly 
affected in this model. In particular, progenitors of bina- 
ries with black holes that start their evolution with stars 
of rather unequal masses are subject to evolving through 
the common envelope and merging. Basically, it is pre- 
dicted that if, in fact, there is no survival in a common 
envelope with Hertzsprung gap donors, then there are al- 
most no BH-BH and BH-NS binaries in our Galaxy within 
the LISA frequency band. In such a case, out of all double 
compact objects, only close NS-NS binaries would have a 
chance to show up in the LISA datastream. 

Although the number of potential LISA NS-NS binaries 
is reduced, it is still quite significant (^ 8 x 10''). The 
physical properties of LISA NS-NS binaries in Model B 
are not much different from Model A, so it is expected 
that the GR signal will have similar shape and will only 
be reduced in strength. For comparison with model A we 
show physical properties of model B potential LISA double 
compact objects in Figure [2j 

3.2. Spatial Distribution 

We present the cumulative distribution (after propaga- 
tion) of distances of DCO from the Earth in Figure [3] 
The distributions corresponding to DCO originating in the 
three components: disk, bulge and halo are shown sepa- 
rately. The distribution of binaries in each component is 
determined by their initial positions and velocities, as well 
as the distribution of merger times. If the mean merger 
time is longer than a few orbits than the distribution of 
binaries may relax to a state determined only by the distri- 
bution of energies and angular momenta of DCOs in each 
component. Since the merger time distribution is typically 
oc the mean merger time is typically half the age of the 
Galaxy, i.e. much longer than a few orbits in the Galactic 
potential . Thus the distribution of DCOs in each compo- 
nent is close to equilibrium in the potential. The energies 
and angular momenta of the binaries are determined by 
their initial positions and velocities. The halo originating 
binaries are born far from the Galactic center and their po- 
tential energies are large, the kick velocities may increase 
them above zero and unbind some of them. There is no 
preferred direction of the angular momentum for these bi- 
naries. Thus their distribution is nearly isotropic around 
the Galactic center. The density profile is almost flat in 
the central 10 — 20kpc, and then is starts to fall off. From 



the Earth we see their number increasing initially as the 
distance cube and then the dependence flattens out beyond 
« 20kpc. The bulge originating binaries have low poten- 
tial energies, and adding the kick velocities unbinds only 
a small fraction of them. Their angular momenta are not 
large since they are all born close to the Galactic center. 
Thus they remain concentrated in the bulge. The cumu- 
lative number as seen from the Earth is determined by 
the small number of escaping binaries at distances smaller 
than « 7kpc, than it goes up very fast when the Galactic 
center is probed. Above « 9kpc the distribution flattens 
out. The binaries originating in the disk present a more 
complicated case. Their potential energies have a wide 
distribution, depending on the distance from the Galactic 
center, however, the disk itself is a potential well. The 
addition of kick velocities increases the potential energy. 
The angular momenta of the disk binaries on their Galac- 
tic orbits are initially aligned perpendicularly to the disk 
plane, yet after the propagation the kick velocities tend 
to partially isotropize them. Thus, they form a disk with 
scale height of a « 0.2kpc and the radius of ~ 15kpc. 
The distribution as seen from the Earth initially probes 
this disk for small distances and the number of DCOs in- 
creases with a cube of distance, yet this dependence is 
only seen for the closest few tens of binaries. For the dis- 
tances from a few hundred pc to ~ lOkpc the number of 
DCOs increases roughly as the distance squared because a 
flat distribution of sources is probed. For larger distances 
(> 15kpc) the dependence flattens out because we reach 
the end of the disk. The disk originating binaries are the 
main contribution to the sources observable from Earth 
for all distances. 

3.3. Gravitational Radiation Signature 

The gravitational wave signature for a 1-year observa- 
tion has been calculated for both models. The population 
has been separated into the different binary types in Fig- 
ure |4] and compared with a realization of the LISA noise. 
In Figure [51 we separate out the different components of 
the Galaxy. From this, we see that in both models, the 
disk binaries are the only ones to rise above the LISA noise. 

We are interested in determining if LISA observations 
can distinguish between these two models. Although we 
do not have an algorithm for identifying eccentric stellar 
mass compact object binaries in the LISA data stream, we 
can still approximate the outcome of such an algorithm. 
It is possible that the coherent signals from the many har- 
monics of eccentric binaries can be summed to boost a 
signal above the LISA noise even though each individual 
harmonic lies buried in the noise (Benacquista 2001, 2002, 
Larson & Hellings 2009). We assume that some sort of 
matched-filter will be used, and so we can describe the de- 
tection statistic as the signal-to-noise ratio (p) defined as 
(Wainstein & Zubakov, 1962): 



HfMf) 



Snif) 



df, 



(20) 



where S'„(/) is the one-sided noise power spectral density. 
For the purposes of calculating the detection statistic, we 
use the discrete analog of Equation [^D] given by Schutz 
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(1997): 

I ~ 2 



and take Sk to be the running average over 1000 bins of the 
power spectrum of the LISA noise given by Equation [19] 
combined with a total Galactic white dwarf binary popula- 
tion generated using Star Track (Ruiter et al. 2009a, b) and 
the DCO population considered in this work. We calcu- 
late hk for each binary in this work using the LISA signal 
simulator described in Section [2731 

If the foreground noise from the confusion-limited signal 
from the Galactic population of close white dwarf bina- 
ries were Gaussian, we could simply assign a reasonable 
threshold value of p as a criterion for detection. However, 
the foreground is not Gaussian (Timpano et al. 2006). 
Racine & Cutler (2007) have studied the impact that the 
non-Gaussianity of the WD- WD foreground has on the de- 
tection of extreme mass ratio inspirals (EMRIs) and super- 
massive black hole inspirals (SMBHs) . Their conclusion is 
that the threshold criterion is essentially unchanged after 
the loudest white dwarf binary signals have been removed 
from the foreground signals. Their result can be explained 
by the fact that the time-scale of the evolution of EMRIs, 
SMBHs, and WD-WDs are all different. Although the 
evolutionary time-scales of stellar mass compact object bi- 
naries are not significantly different from those for white 
dwarf binaries, white dwarf binaries are almost exclusively 
circularized in the LISA band and nearly all compact ob- 
ject binaries analyzed in this work are eccentric. Thus, 
their waveforms differ significantly in the time domain. 
This can be used to distinguish between the two popula- 
tions, and we assume that the results of Racine and Cutler 
(2007) can be applied here. 

Finally, we note that the imposition of a threshold 
signal-to-noise ratio calculated according to Equation [^1] 
can result in detection of binaries for which no individ- 
ual harmonic has a signal-to- noise ratio above 1. Matched 
filtering algorithms have already been demonstrated that 
can successfully detect SMBH signals within mock LISA 
data (Babak et al. 2008) even though the signal never 
exceeds the instrument and white dwarf binary confusion 
noise at any frequency (see Fig 1 of Arnaud et al. 2007). 

Consequently, we impose a threshold of p > 10 as the 
criterion for detection. Using this estimate, we find that 
there are 11 resolvable double compact object binaries in 
model A and only 2 in model B. From Figure [H we see 
that the bulk of the potentially detectable very low fre- 
quency sources in model A are BH-BH binaries and that 
these are dramatically suppressed in model B. This is also 
borne out by the estimates of resolvable binaries where 
the number of resolvable BH-BH binaries drops from 5 in 
model A to zero in model B. Thus, detection of high-mass, 
low-frequency binaries can be used to distinguish between 
these two models. In model A there is a relatively large 
number of BH-BH binaries at low frequencies (at 10"** Hz 
there is a comparable number of NS-NS and BH-BH sys- 
tems; see Fig. [1]). Since at these frequencies binaries have 
relatively large orbital separations, neutron stars that are 
of much lower mass than black holes do not generate a de- 
tectable GR signal, while massive black hole systems do. 
For lower frequencies (^ 10^^ Hz) corresponding to yet 



larger orbital separations, even the most massive BH-BH 
binaries do not generate detectable signal, although they 
dominate over other double compact objects. The proper- 
ties of the potentially resolvable binaries are presented in 
Table [5] and shown in Figures [T] and [5] 

4. COMPARISON WITH PREVIOUS WORK 

The number of detectable compact object binaries that 
we found in this analysis is somewhat smaller than the 42 
detectable systems predicted by Nelemans et al. (2001). 
However, this can be reasonably explained by the differ- 
ences in the overall population numbers between the two 
simulations. The total number of compact object binaries 
predicted in our simulation ranges between 1.0 x 10^ for 
model B and 1.4 x 10® for model A, while Nelemans et 
al. predict 4.0 x 10® such systems. Thus, our prediction 
of ~ 11 systems is comparable with Nelemans et al. It 
should be noted, however, that our estimate is conserva- 
tive as we have assumed no prior removal of any WD- WD 
signals from the calculated value of Sk-, and that we have 
assumed an observation time of 1 year. In reality, we ex- 
pect some removal of the loud WD- WD signals and the 
nominal mission lifetime of LISA is 3 years. On the other 
hand, we have not included any penalties incurred by the 
larger parameter space needed for the analysis of eccentric 
systems. 

5. CONCLUSIONS 

We have calculated, in a self-consistent way, the gravita- 
tional radiation signal of Galactic double compact object 
binaries. The realistic model for Galactic stars was em- 
ployed, then massive stars were followed in their evolution 
to form neutron stars and black holes, some of which were 
formed in binaries. The current position of the binaries 
with neutron stars and black holes relative to Earth was 
obtained and integrated gravitational radiation signal was 
predicted for the LISA detector. 

It was found that several double compact object binaries 
can be detected and resolved in a 1-year observation with 
an instrument like LISA. In model A, we find comparable 
numbers of BH-BH and NS-NS binaries, while no BH-BH 
binaries are observed in model B. This prediction, if true, 
would put a new light on double compact object studies as 
at present we have only detected 9 NS-NS systems. Double 
compact objects with a black hole are yet to be discovered. 
LISA could potentially detect several of these binaries. 
Starting with NS-NS systems, these potential detections 
could answer a number of questions. Do nearly all NS-NS 
systems host low mass neutrons stars (~ 1.35 Mq) as it 
seems to appear in a known sample discovered through 
radio-pulsar surveys or is it only some observational bias? 
Galactic merger rates seem to appear consistent with the 
notion that NS-NS are short GRB progenitors. But would 
this be confirmed by a broader sample detected by LISA? 
Detection of any double compact objects with black holes 
would be a feat in and of itself. Although a number of these 
systems are expected from evolutionary calculations, these 
are usually burdened with rather large uncertainties. De- 
tections of BH-BH binaries would provide very valuable 
tests for evolutionary codes that are used in number of 
studies. Knowledge of black hole numbers and their basic 
properties can put direct constraints on evolution of the 
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most massive stars, their progenitors. In particular, one 
such example was presented in our study. Depending on 
a common envelope model the number of BH-BH systems 
that could be detected with LISA vanishes. So just by 
their presence (or lack of thereof) in the signal one could 
test an evolutionary phase that is virtually undetectable 
due to its shortness (~ 10"^ — 10'' yr). And it it interesting 
to note, that although the common envelope was proposed 
more than 30 years ago (Paczynski 1976), there is still a 
surprising lack of understanding of the outcome — and the 
outcome is crucial for formation of Type la progenitors 
(close white dwarf binaries) or any type of X-ray binaries 
that are currently observed even in external galaxies by 
Chandra or XMM, just to name a few examples. 

Although, double compact objects will be most proba- 
bly discovered first in gravitational radiation with ground- 
based interferometers (e.g., LIGO, VIRGO) they will be 
detected only in the inspiral and merger phases, prevent- 
ing a measurement of their orbital parameters (e.g., sepa- 
rations and eccentricities). Additionally, the merger rates 
are too small to expect detections in the Galaxy, and the 



observations will be made for a wide spectrum of local en- 
vironments in which mergers are taking place, making it 
harder to connect the result to a particular evolutionary 
model that by nature depends on local conditions (e.g., 
metallicity, IMF or age of population). In this context, 
LISA may provide a potentially very unique sample of 
double compact object binaries in the most studied en- 
vironment of our own Galaxy and in the same providing 
the most robust constraints on evolution of massive stars. 
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Table 1 

Number of Double Compact Binaries in the Galaxy*" 



NS-NS BH-NS BH-BH 



ALL: disk 


281091 - 


238245 


45075 - 


17729 


621812 - 


424501 


ALL: bulge 


58658 - 


52520 


10347 - 


- 4388 


147260 - 


106170 


ALL: halo 


6351 


6342 


2845 - 


2776 


37681 - 


36373 


ALL: total 


346100 - 


297107 


58267 - 


24893 


806753 - 


567044 


LISA: disk 


92253 - 


58633 


4029 


371 


6903 


- 6 


LISA: bulge 


12853 


8700 


382 - 


- 48 


417- 


- 


LISA: halo 


298 - 


291 


50 - 


39 


36 - 


25 


LISA: total 


105404 - 


- 67624 


4461 - 


- 458 


7356- 


- 31 



''Numbers are given for two models: Model A - Model B. ALL - 
all binaries in a given category. LISA - binaries with orbital periods 
and eccentricities that are likely to contribute to GR signal in LISA 
frequency band. 



Table 2 

Properties of resolvable Double Compact Binaries in the 

Galaxy 



Model 


Type 


-^chirp 

(Mo) 


/orb 

(mHz) 


e 


P 


A 


BH-BH 


5.35 


0.284 


0.038 


22.7 


A 


BH-BH 


5.77 


0.455 


0.005 


22.1 


A 


BH-BH 


4.01 


0.448 


0.011 


17.3 


A 


BH-BH 


7.69 


0.097 


0.006 


15.2 


A 


BH-BH 


7.46 


0.279 


0.006 


11.5 


A 


NS-NS 


1.16 


1.449 


0.167 


53.4 


A 


NS-NS 


1.24 


0.187 


0.884 


32.1 


A 


NS-NS 


1.16 


0.760 


0.346 


17.3 


A 


NS-NS 


1.15 


0.956 


0.044 


13.2 


A 


NS-NS 


1.16 


0.843 


0.106 


12.2 


A 


NS-NS 


1.15 


0.990 


0.069 


10.0 


B 


NS-NS 


1.16 


0.760 


0.346 


17.3 


B 


NS-NS 


1.15 


0.990 


0.069 


10.0 



10 
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Fig. 1. — Potential LISA Galactic population of double compact objects (Model A). The top panel shows the distribution 
of orbital frequency, middle panel the distribution of eccentricities, while the bottom panel has the distribution of chirp 
masses. Note the dominance of NS-NS binaries. At the top of each panel we show Galactic resolved double compact 
objects. Since this population is very small (6 NS-NS and 5 BH-BH binaries) we show all individual resolved sources 
(details in Tabled]). 
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Fig. 2. — Same as Fig. [T]but for model B. The NS-NS dominance is even more pronounced in this model. Note that there 
are only 2 resolved sources and they are both NS-NS binaries. 
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Fig. 3. — The cumulative distribution of the current distances of simulated binaries from the Earth (model A). The 
distances are obtained after propagation in the Galactic potential; the binaries have spread out from their birth sites 
due to natal kicks. We show the three Galactic components separately. Note that the Galactic center is at a distance of 
~ 8kpc from Earth. 
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Fig. 4. — Strain spectral densities of the double compact object binaries separated by binary type and compared with 
the full spectrum containing the double compact objects, the WD- WD signal and the LISA noise. The line is a running 
average over the WD- WD and LISA noise spectrum. The left panel shows model A and the right panel shows model B. 
Note the strong suppression of the BH-BH binaries in model B, and the absence of any detectable signals below ^ 0.4 
mHz. 
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Fig. 5. — Strain spectral densities of the double compact object binaries separated by Galaxy component and compared 
with the full spectrum containing the double compact objects, the WD- WD signal and the LISA noise. The line is a 
running average over the WD- WD and LISA noise spectrum. The left panel shows model A and the right panel shows 
model B. Note the dominance of the disk binaries in both models, but the significant absence of low frequency binaries in 
model B. 



